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ABSTRACT A 32-nucleotide (nt) RNA motif located at the 3’ end of the transmissible gastroenteritis coronavirus (TGEV) ge- 
nome was found to specifically interact with the host proteins glutamyl-prolyl-tRNA synthetase (EPRS) and arginyl-tRNA syn- 
thetase (RRS). This RNA motif has high homology in sequence and secondary structure with the gamma interferon-activated 
inhibitor of translation (GAIT) element, which is located at the 3’ end of several mRNAs encoding proinflammatory proteins. 
The GAIT element is involved in the translation silencing of these mRNAs through its interaction with the GAIT complex (EPRS, 
heterogeneous nuclear ribonucleoprotein Q, ribosomal protein L13a, and glyceraldehyde 3-phosphate dehydrogenase) to favor 
the resolution of inflammation. Interestingly, we showed that the viral RNA motif bound the GAIT complex and inhibited the in 
vitro translation of a chimeric mRNA containing this RNA motif. To our knowledge, this is the first GAIT-like motif described in 
a positive RNA virus. To test the functional role of the GAIT-like RNA motif during TGEV infection, a recombinant coronavirus 
harboring mutations in this motif was engineered and characterized. Mutations of the GAIT-like RNA motif did not affect virus 
growth in cell cultures. However, an exacerbated innate immune response, mediated by the melanoma differentiation-associated 
gene 5 (MDA5) pathway, was observed in cells infected with the mutant virus compared with the response observed in cells in- 
fected with the parental virus. Furthermore, the mutant virus was more sensitive to beta interferon than the parental virus. All 
together, these data strongly suggested that the viral GAIT-like RNA motif modulates the host innate immune response. 


IMPORTANCE The innate immune response is the first line of antiviral defense that culminates with the synthesis of interferon 
and proinflammatory cytokines to limit virus replication. Coronaviruses encode several proteins that interfere with the innate 
immune response at different levels, but to date, no viral RNA counteracting antiviral response has been described. In this work, 
we have characterized a 32-nt RNA motif located at the 3' end of the TGEV genome that specifically interacted with EPRS and 
RRS. This RNA motif presented high homology with the GAIT element, involved in the modulation of the inflammatory re- 
sponse. Moreover, the disruption of the viral GAIT-like RNA motif led to an exacerbated innate immune response triggered by 
MDAS, indicating that the GAIT-like RNA motif counteracts the host innate immune response. These novel findings may be of 
relevance for other coronaviruses and could serve as the basis for the development of novel antiviral strategies. 
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“oronaviruses (CoVs) are enveloped, single-stranded, (5) and the recently emerged (April 2012) Middle East respiratory 


“= positive-sense RNA viruses that belong to the Coronaviridae 
family within the Nidovirales order (1). CoVs are vertebrate 
pathogens responsible mainly for respiratory and enteric infec- 
tions in a wide range of animals and human (2). Among the high 
diversity of CoVs infecting animal species, transmissible gastroen- 
teritis virus (TGEV) is of special relevance in pigs, causing a life- 
threatening disease with important economic losses (3). In hu- 
mans, CoV infections have been historically associated with mild 
upper respiratory tract diseases (4). However, the identification of 
the severe acute respiratory syndrome CoV (SARS-CoV) in 2003 
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syndrome CoV (MERS-CoV) (6), both causing acute pneumonia 
and even death, redefined historic perceptions and potentiated the 
relevance of CoVs as important human pathogens. Therefore, un- 
derstanding the molecular basis of CoV replication and pathogen- 
esis will allow the development of effective strategies to prevent 
and control CoV infections. 

CoVs contain the largest known RNA genome among RNA 
viruses, consisting of a plus-sense, 5'-capped, and polyadenylated 
RNA molecule of 27 to 31 kb in length (2, 7). The first two-thirds 
of the genome carries the replicase gene, which comprises two 
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overlapping open reading frames (ORFs) named 1a and 1b. Both 
ORFs are directly translated from the viral genome, resulting in 
two large polyproteins (ppla and pplab), which are autoproteo- 
litically cleaved to release the replication-transcription complex 
components (8). The 3’ one-third of the genome includes the 
genes encoding the structural spike (S), envelope (E), mem- 
brane (M), and nucleocapsid (N) proteins, as well as the genus- 
specific proteins characteristic of each CoV, which are ex- 
pressed from a nested set of 3’-coterminal subgenomic mRNAs 
(sgmRNAs) (2, 7, 8). 

CoV replication and transcription are complex processes that 
require the specific recognition of RNA cis-acting elements lo- 
cated at the ends of the viral genome (9) and are mediated by a 
huge protein complex encoded by the viral replicase gene together 
with the participation of cellular proteins (2, 7, 9). To date, several 
cellular proteins implicated in CoV replication have been identi- 
fied on the basis of their interaction with replicase proteins (10, 
11) or with cis-acting elements required for viral genome replica- 
tion (12, 13). Recently, we identified 10 cellular proteins interact- 
ing with the TGEV genome ends by RNA affinity chromatography 
and mass spectrometry analysis (13). These proteins included the 
polypyrimidine tract-binding protein (PTB), which preferentially 
interacted with the 5’ end of the genome, and a set of 9 proteins, 
including several heterogeneous nuclear ribonucleoproteins (hn- 
RNPs) (Al, A2B1, AO, Q, and VU), the translation factors glutamyl- 
prolyl-tRNA synthetase (EPRS), arginyl-tRNA synthetase (RRS), 
and poly(A)-binding protein (PABP), and the p100 transcrip- 
tional coactivator, which bound preferentially to the 3’ end of the 
genome. Silencing studies demonstrated a negative role in viral 
RNA synthesis for PTB (9) and a positive role for the 3’-end- 
interacting proteins hnRNP Q, EPRS, RRS, and PABP (13) (S. 
Marquez-Jurado, L. Enjuanes, and F. Almazan, unpublished 
data). 

The innate immune response is the first line of antiviral de- 
fense. During infection, viral components are detected through a 
set of pathogen recognition receptors (PRRs) (14, 15), including 
the Toll-like receptors (TLRs), the retinoic acid-inducible gene I 
product (RIG-I), and the melanoma differentiation-associated 
gene 5 (MDAS) protein. This recognition triggers a signaling cas- 
cade that culminates with the synthesis of interferon (IFN) and 
proinflammatory cytokines to limit viral replication and spread 
(16). For RNA viruses, double-stranded RNAs (dsRNAs), which 
are present in the viral genome as stem-loop structures or are 
generated during viral replication, are the most important 
pathogen-associated molecular pattern (PAMP) recognized by 
cellular PRRs (15, 17). In the case of CoVs, MDAS has been de- 
scribed as the main cytoplasmic PRR sensing the viral RNAs (18). 

To evade the innate immune response, viruses have developed 
an extensive and highly diverse range of countermeasures to im- 
pair the antiviral response, including blocking different compo- 
nents of the innate immune pathways or hiding dsRNAs to make 
them inaccessible to sensors (14). In this sense, CoVs encode sev- 
eral proteins that interfere with type-I IFN and proinflammatory 
cytokine production at different levels (15). Similar to that of 
other positive-stranded RNA viruses, CoV replication takes place 
in the cytoplasm, leading to the generation of dsRNA molecules 
that trigger the host antiviral innate immune response. To over- 
come this response, CoVs have developed several strategies, in- 
cluding the induction of double-membrane vesicles (DMVs) (14), 
proposed to hide viral dsRNAs from innate sensors during ge- 
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nome replication. In addition to these dsRNA molecules, most 
likely, the viral genome and the ssmRNAs of CoVs could contain 
highly structured RNA domains that could be recognized by the 
cell guard sensors (19), as described for other viruses, like hepatitis 
C virus (20), reovirus (21), coxsackievirus (22), West Nile virus, 
and yellow fever virus (23). However, the nature of these highly 
structured RNA domains and the molecular mechanisms used by 
CoVs to avoid their recognition by the cellular sensors are still 
unknown. In addition, other viruses, such as adenoviruses, tran- 
scribe two adenovirus-associated RNA I (VAI) which fold into 
stem-loop structures that are involved in the blocking of some 
components of the innate immune pathway, such as the protein 
kinase R (PKR) (24). Whether RNA viruses, such as CoVs, contain 
similar RNAs interfering with the innate immune response re- 
mains to be determined. 

In this work, a novel RNA domain located at the 3’ end of the 
TGEV genome that specifically interacted with the aminoacyl- 
tRNA synthetases EPRS and RRS has been identified and func- 
tionally characterized. This RNA motif presented high homology 
in sequence and secondary structure with the gamma interferon- 
activated inhibitor of translation (GAIT) element that is located at 
the 3’ end of several mRNAs encoding proinflammatory proteins. 
This element is involved in the translation silencing of these mR- 
NAs to favor the resolution of inflammation. The functional role 
of the viral GAIT-like RNA motif during TGEV infection was 
further analyzed by mutational analysis using a TGEV infectious 
clone, demonstrating that the viral GAIT-like RNA motif influ- 
enced the innate immune response by affecting the MDA5 path- 
way. To our knowledge, this is the first time that a CoV RNA 
domain has been directly involved in the modulation of the innate 
immune response. 


RESULTS 
Relevance of the poly(A) tail in the binding of cellular proteins 
to the 3’ end of the TGEV genome. In a previous study, nine 
cellular proteins were found to bind preferentially to the 3’ end of 
the TGEV genome, including several hnRNPs (A1, A2B1, AO, Q, 
and U), the translation factors EPRS, RRS, and PABP, and the 
p100 transcriptional coactivator (13). The specific RNA domain 
interacting with these cellular proteins was further analyzed. In a 
first approach, the role of the poly(A) tail in the binding of these 
proteins with the 3’ end of the genome was studied. To this end, 
three RNA constructs, one containing the first 504 nucleotides 
(nt) of the 5’ end of the genome and the others containing the last 
493 nt of the 3’ end with or without the poly(A) tail, were gener- 
ated by in vitro transcription, labeled with biotin, and used as baits 
for RNA affinity protein isolation (Fig. 1). The RNA baits were 
incubated with cytoplasmic cell extracts from infected human 
liver-derived Huh-7 cells, and the pulled-down proteins were de- 
tected by Western blotting by using specific antibodies. Huh-7 
cells were used instead of swine testis (ST) cells, because Huh-7 
cells are susceptible to TGEV and because of the availability of 
validated antibodies against human proteins (13). As expected, 
the PABP specifically bound to the 3’ end containing the poly(A) 
tail (Fig. 1). In contrast, all the remaining proteins interacted with 
the 3’ end both with and without the poly(A) tail (Fig. 1), indicat- 
ing that EPRS, RRS, p100, and the hnRNPs Q, Al, U, and AO did 
not require the poly(A) tail to interact with the genome 3’ end. 
Mapping of TGEV 3’-end RNA domains interacting with cel- 
lular proteins. To map the 3’-end RNA domains interacting with 
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FIG 1 Relevance of the poly(A) tail in the binding of the 3’-end-interacting 
cellular proteins. Proteins (600 wg) from cytoplasmic extracts of TGEV- 
infected Huh-7 cells were pulled down without RNA (—) or with TGEV ge- 
nome ends 5’ (5’) and 3’ without (3) and with (3’ A) poly(A) tail, separated in 
SDS-PAGE gels, and analyzed by immunoblotting with specific antibodies 
against PABP, EPRS, RRS, p100, and the hnRNPs (Q, Al, U, and AO). 


= — hnRNP AO 


EPRS, RRS, p100, and the hnRNPs Q, Al, U, and AO, the last 
493 nt of the TGEV genome were divided into six overlapping 
RNA fragments (F1 to F6) that were used as baits in pulldown 
assays (Fig. 2A). These RNA fragments were designed taking into 
consideration the predicted secondary structure of the 3’ end of 
the TGEV genome and the information available about the 3'-end 
cis-acting elements involved in CoV replication. Fragments F1 and 
F2, covering nt 28087 to 28207 and 28147 to 28277, respectively, 
contained four predicted stem-loops for which no function has 
been associated. Fragments F3, F4, and F5, spanning nt 28269 to 
28388, 28358 to 28481, and 28289 to 28426, respectively, covered 
a region containing RNA structures similar to the bulged stem- 
loop (BSL) and the pseudoknot (PK), which have been described 
to be important for genome replication in mouse hepatitis virus 
(MHV) and bovine CoV (BCoV) (25-27). Finally, fragment F6, 
covering nt 28469 to 28580, contained the hypervariable region 
(HVR) and the highly conserved octanucleotide (OCT) (5'’ GGA 
AGA GC 3’), originally described for MHV (27). The biotin- 
labeled RNA fragments, generated by in vitro transcription, were 
used as baits in pulldown assays, together with the full-length 5’ 
and 3’ ends of the genome as specificity controls. In addition to the 
expected preferential binding of all proteins to the 3’ end of the 
genome, a specific binding to F2 was observed for EPRS and RRS, 
and a preferential binding to F4 was found for hnRNP Al 
(Fig. 2B). In the case of hnRNPs Q and AO, these proteins were 
preferentially pulled down with F2 and F6 and with F2 and F4, 
respectively (Fig. 2B). In contrast, no preferential binding for any 
fragment was detected in the case of hnRNP U and p100. Due to 
the specific interaction of EPRS and RRS with F2, and the positive 
role of these proteins in CoV RNA synthesis (13) (Marquez- 
Jurado et al., unpublished), this RNA fragment was selected to 
further delimitate the RNA domains interacting with both 
aminoacyl-tRNA synthetases. 

Identification of the 3’-end RNA domains interacting with 
EPRS and RRS. To delimitate the RNA domains within F2 inter- 
acting with EPRS and RRS, this fragment was divided into three 
smaller RNA fragments (F2.1, F2.2, and F2.3) according to its 
predicted secondary structure (Fig. 3A). The stability of the pre- 
dicted harping loop structures present in these RNA fragments 
was equivalent to that of the corresponding harping loop struc- 
tures presented in fragment F2 and in the full-length 3’ end of the 
genome. These fragments were chemically synthesized, labeled 
with biotin, and used as baits in pulldown assays, together with F2 
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and F3 as specificity controls. In addition to the expected binding 
to F2, a preferential binding to F2.2 was observed for both EPRS 
and RRS (Fig. 3A). 

The F2.2 RNA fragment was subsequently divided into four 
new fragments (F2.2L, F2.2R, F2.2U, and F2.2D), according to its 
sequence and predicted secondary structure that consists of a 
complex stem-loop with two internal bulges (Fig. 3B). F2.2L and 
F2.2R were comprised of the 5’ and 3’ halves of the fragment F2.2 
sequence, respectively, while F2.2U contained the upper stem of 
the structure and F2.2D the two internal bulges. In addition to the 
expected binding to F2.2, both aminoacyl-tRNA synthetases pref- 
erentially bound to fragment F2.2L, which constitutes a 32-nt 
RNA domain located 410 nt upstream of the 3’ end of the viral 
genome, within the coding sequence of gene 7 (Fig. 3C). The se- 
quence of this RNA domain showed high identity with the GAIT 
element, present in the 3’ untranslated region (UTR) of several 
mRNAs coding for proinflammatory proteins, such as the cerulo- 
plasmin (Cp) (Fig. 3C) (28, 29). It has been described that IFN-y 
induces the assembly of the GAIT complex, constituted by EPRS, 
hnRNP Q, the ribosomal protein L13a, and glyceraldehyde 
3-phosphate dehydrogenase (GAPDH). This protein complex in- 
teracts with the GAIT element present in different mRNAs encod- 
ing proinflammatory proteins, inhibiting their translation and 
therefore the harmful accumulation of inflammatory products 
(29, 30). Interestingly, the secondary structure of the GAIT ele- 
ment, consisting of a stem-loop with an internal asymmetric 
bulge, was also predicted for the viral RNA domain, including the 
A and U residues located in the asymmetric internal bulge, which 
have been reported to be critical for the function of the cellular 
GAIT element (Fig. 3C) (28, 29). Overall, these data suggest that 
the viral RNA domain, expanding fragment F2.2L, could consti- 
tute a GAIT-like RNA motif. 

Analysis of the interaction of EPRS and RRS with viral RNAs 
during infection. The interaction of EPRS and RRS with viral 
RNAs in the context of TGEV infection was analyzed by RNA 
immunoprecipitation assays. To this end, cytoplasmic extracts of 
TGEV-infected ST cells were immunoprecipitated with EPRS- 
and RRS-specific antibodies, and the presence of the viral gnome 
and sgmRNAs in the immunoprecipitated RNA-protein com- 
plexes was analyzed by reverse transcription-quantitative PCR 
(RT-qPCR) by using specific oligonucleotides. In contrast to the 
samples immunoprecipitated with the negative-control anti- 
green fluorescent protein (GFP) antibody, the viral genomic RNA 
and sgmRNAs were specifically detected in the RNA-protein com- 
plexes immunoprecipitated with anti-EPRS or anti-RRS antibod- 
ies (Fig. 4). These data indicated that both aminoacyl-tRNA syn- 
thetases interacted with the viral genome and the sgmRNAs 
during infection, as expected from the presence of the GAIT-like 
RNA motif in these viral RNAs. 

In vitro analysis of the interaction of EPRS and RRS with the 
viral GAIT-like RNA motif. The pulldown assays did not discrim- 
inate whether the interaction of EPRS and RRS with the viral 
GAIT-like RNA motif was direct or mediated by other cellular or 
viral proteins. To clarify this issue, this interaction was analyzed by 
electrophoretic mobility shift assay (EMSA) by using purified pro- 
teins and the biotin-labeled GAIT-like RNA motif (Fig. 5). Puri- 
fied GAPDH, a protein of the GAIT complex that does not directly 
interact with the GAIT element (30), and the DDK/Myc tag pres- 
ent in the purified EPRS were used as negative controls. A specific 
interaction of the TGEV GAIT-like RNA motif with the purified 
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FIG 2 Mapping of the RNA domains interacting with the 3’-end-interacting cellular proteins. (A) Predicted secondary structure of the 493 nt from the TGEV 
genome 3’ end, obtained using the Mfold 3.4 software (49). The sequence covered by the six overlapping fragments (F1 to F6) and the higher-order RNA 
structures similar to the bulge stem-loop (BSL), the pseudoknot (PK), the hypervariable region (HVR), and the conserved octanucleotide (OCT) are indicated. 
(B) Identification of the 3’-end domains interacting with the selected cellular proteins. Proteins (600 jg) from cytoplasmic extracts of TGEV-infected Huh-7 cells 
were pulled down with the six overlapping fragments (1 to 6) covering the last 493 nt of the TGEV genome or with TGEV genome ends 5’ and 3’ without poly(A) 
or without RNA (—) as specificity controls. The pulled down proteins were separated in SDS-PAGE gels together with 50 jg of the initial protein extract (I) and 
analyzed by immunoblotting with specific antibodies against EPRS, RRS, hnRNP Q, hnRNP Al, hnRNP AO, hnRNP U, and p100. To simplify the information 
of the figure, an empty lane was removed in the hnRNP AO panel, in which the dashed line indicates the splicing site. The molecular masses in kDa of the analyzed 


proteins are indicated at the left. 


EPRS but not with the DDK/Myc tag or GADPH was observed 
(Fig. 5). Unfortunately, in the case of RRS, the purified protein did 
not enter into the gels used for the EMSA, and nonconclusive data 
were obtained using this technique (data not shown). These data 
indicated that, at the least, the EPRS interacted directly with the 
TGEV GAIT-like RNA motif in vitro, similar to what has been 
described for the cellular GAIT element (29, 31). 

Invitro functional analysis of the viral GAIT-like RNA motif. 
Based on the functional role of the cellular GAIT element on 
translation silencing (29) and the direct binding of EPRS to the 
viral GAIT-like RNA motif, the potential inhibitory effect of this 
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viral RNA motif on translation was analyzed in vitro. To this end, 
three chimeric mRNAs expressing the luciferase protein, contain- 
ing either the Cp GAIT element, the TGEV GAIT-like RNA motif, 
or no GAIT element in the 3’ UTR, were in vitro transcribed, and 
their translation efficiency was analyzed in the presence or in the 
absence of the GAIT complex by using a rabbit reticulocyte lysate 
system (Fig. 6). It has been described that the GAIT complex is 
assembled in monocytic cells after treatment with IFN-y (30). 
Therefore, cytosolic extracts of human U937 monocytic cells un- 
treated or treated with IFN-y were used in the in vitro translation 
reactions as a source of the GAIT complex. As expected, transla- 
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FIG3 Mapping of the RNA domains interacting with EPRS and RRS. (A) Mapping of the region within F2 interacting with ERPS and RRS. F2 was divided into 
three fragments (F2.1, F2.2, and F2.3) that were used as bait in pulldown assays, together with fragments F2 and F3 or without RNA (—) as specificity controls. 
The presence of EPRS and RRS in the initial extract (I) and in the pulled-down fractions was analyzed by Western blotting. To simplify the information in the 
figure, an empty lane was removed in the EPRS panel, in which the dashed line indicates the splicing site. The molecular masses in kDa of EPRS and RRS are 
indicated at the left. (B) Identification of the F2.2 domain interacting with EPRS and RRS. F2.2 was divided into four fragments (F2.2L, F2.2R, F2.2U, F2.2D), and 
their interactions with EPRS and RRS were analyzed by pulldown assays and Western blotting. The same experiment was performed with fragment F2.2 or 
without RNA (—) as specificity controls. The molecular masses in kDa of EPRS and RRS are indicated at the left. (C) Sequence analysis of fragment F2.2L. The 
bar on the top represents the TGEV genome, in which the different genes (ORF la, ORF 1b, S, 3a, 3b, E, M, N, and 7 genes), the leader sequence (L), and the 3’ 
UTR are illustrated as boxes. The sequence alignment of the viral fragment F2.2L with the GAIT element of the ceruloplasmin (Cp GAIT) and their respective 
secondary structures, in which the A and U residues critical for the function of the GAIT element are outlined in green, are shown. Numbers above the sequence 
alignment and in the secondary structure indicate the genome position of fragment F2.2L. 


tion of the luciferase gene was detected in all cases in the absence of 
the GAIT complex (Fig. 6). In contrast, a significant inhibition of 
luciferase translation was observed from the mRNA constructs 
containing either the cellular (75%) or the viral (50%) GAIT ele- 
ment in the presence of the GAIT complex (Fig. 6). Interestingly, 
when increasing amounts of either the cellular or the viral GAIT 
elements were added to the translation reaction mixtures as com- 
petitors, a reversion of the luciferase translation inhibition was 
observed (data not shown), indicating that the translation inhibi- 
tion is due to the presence of the GAIT element. These data indi- 
cated that, similar to the cellular GAIT element, the viral GAIT- 
like RNA motif had an inhibitory effect on in vitro translation, 
probably by interaction with the GAIT complex. 

Generation and characterization of a recombinant TGEV 
harboring mutations in the GAIT-like RNA motif. Once the in- 
hibitory effect on translation of the viral GAIT-like RNA motif 
was determined in vitro, the functional role of this RNA motif was 
studied by mutational analysis in the context of the TGEV infec- 
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FIG4 Analysis of the interaction of EPRS and RRS with viral RNAs by RNA 
immunoprecipitation assay. Cytoplasmic extracts of TGEV-infected ST cells 
were immunoprecipitated with specific antibodies for EPRS (a@-EPRS), RRS 
(a@-RRS), and GFP (a-GFP). The presence of viral genomic RNA (gRNA) and 
sgmRNAs (S, 3a, M, E, N, and 7 genes) in the immunoprecipitated RNA- 
protein complexes was quantified by RT-qPCR and expressed as the fold 
change of the RNA levels relative to the RNA levels in the samples immuno- 
precipitated with GFP antibody. Error bars represent the standard deviations 
from three independent experiments. 
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FIG5 Study of the interaction between EPRS and the TGEV GAIT-like RNA 
motif by RNA EMSA. Purified EPRS, DDK/Myc tag, and GAPDH proteins 
were incubated with the 3’-biotinylated GAIT-like RNA motif riboprobe 
(GAIT-like). The RNA-protein complexes were resolved by electrophoresis on 
a nondenaturing 4% polyacrylamide gel, and the biotinylated riboprobe was 
detected using the BrightStar BioDetect kit (Ambion). The arrow indicates the 
shifted protein-RNA complex. 


tion. Based on previously reported functional data of the Cp GAIT 
element (28, 29), a TGEV recombinant virus harboring mutations 
in the GAIT-like RNA motif (rTGEV-mutG) was engineered by 
reverse genetics using a TGEV infectious clone (Fig. 7A) (32). The 
rTGEV-mutG virus contained the nucleotide substitutions A to G 
(genomic position 28178), U to C (genomic position 28181), and 
C to U (genomic position 28190), which were predicted to par- 
tially disrupt the stem structure and the internal asymmetric 
bulge. As the GAIT-like RNA motif is located within the coding 
region of gene 7, the mutations introduced were silent and main- 
tained the amino acid sequence of protein 7. The rTGEV-mutG 
virus was successfully rescued and the mutations introduced were 
present in the recovered mutant virus, even after 20 serial passages 
in tissue culture of a plaque-purified virus, indicating that they 
were stably maintained in the rTGEV-mutG genome. In addition, 
no significant differences in the plaque size and the growth kinet- 
ics, at both low (0.05) and high (3) multiplicities of infection 
(MOI), were observed between the mutant and wild-type (WT) 
viruses (Fig. 7B), suggesting that the viral GAIT-like RNA motif 
was not essential for virus replication in cell cultures. 

The generation of rTGEV-mutG was based on the assumption 
that the introduced mutations abolish the binding of EPRS and 
RRS to the viral-GAIT-like RNA motif and, therefore, inactivate 
the putative function of this RNA motif. To confirm that this was 
the case, the binding of both aminoacyl-tRNA synthetases to ei- 
ther the native or the mutated GAIT-like RNA motif (see Table 2, 
F2.2L or F2.2L-mutG, respectively) was analyzed using pulldown 
assays (Fig. 8). As expected, the mutated viral GAIT-like RNA 
motif was not able to pull down both aminoacyl-tRNA syntheta- 
ses, indicating that EPRS and RRS did not bind to the mutated 
viral GAIT-like RNA motif (Fig. 8). 

All together, these data indicated that the viral GAIT-like RNA 
motif was dispensable for virus replication in cell culture. 

Effect of GAIT-like mutations on the innate immune re- 
sponse. Due to the dispensability of the GAIT-like RNA motif for 
virus replication in cell culture, and considering the reported role 
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FIG 6 In vitro analysis of the inhibitory effect on translation of the TGEV 
GAIT-like RNA motif. Three chimeric luciferase reporter mRNA constructs, 
containing at the 3’ UTR no GAIT element (No GAIT), the Cp GAIT element 
(Cp GAIT), or the TGEV GAIT-like RNA motif (TGEV GAIT), were subjected 
to in vitro translation in a rabbit reticulocyte lysate in the absence (blue bars) or 
in the presence (red bars) of the GAIT complex. The translation of these chi- 
meric constructs was quantitated by chemiluminescence and expressed as the 
percentage of translation of each construct in the presence of the GAIT com- 
plex compared to that in the absence of the GAIT complex. Error bars repre- 
sent the standard deviations from three independent experiments. ***, P < 
0.001. 


of the Cp GAIT element modulating the inflammatory response 
(30, 33), the potential role of the TGEV GAIT-like RNA motif in 
the modulation of the host inflammatory response was further 
studied. To this end, ST cells were infected at an MOI of 3 with the 
rTGEV-WT and rTGEV-mutG viruses, and the mRNA levels of a 
set of proinflammatory cytokines (tumor necrosis factor [TNF], 
CCL2, and CXCL9) were analyzed at different hours postinfection 
(hpi) by RT-qPCR using specific TaqMan assays. As a control of 
viral replication, the levels of genomic RNA and sgmRNA 7 were 
analyzed. In addition, to discard a general transcription up- or 
downregulation during virus infection that could affect the ex- 
pression of all cellular genes, the mRNA levels of the transforming 
growth factor B (TGF-) gene, which is involved in the immune 
response, were also evaluated as a control. As expected, all the 
analyzed genes increased their expression during the WT virus 
infection (Fig. 9). Interestingly, all of them were significantly up- 
regulated (from 2- to 20-fold, depending on the gene) in cells 
infected with the rTGEV-mutG virus compared with in cells in- 
fected with the WT virus (Fig. 9). In contrast, no significant dif- 
ferences were detected in the levels of viral genomic RNA, sgm- 
RNA 7, and TGF-B mRNA in cells infected with either the WT or 
mutant viruses (Fig. 9), discarding that proinflammatory RNA 
upregulation was due to a difference in virus replication or to a 
general transcription upregulation during r[GEV-mutG infec- 
tion. 

The host innate immune response led to both proinflamma- 
tory cytokines and type I IFN production. Therefore, IFN-B and 
IFN-stimulated genes (ISGs) (MDA5 and RIG-I) production was 
also analyzed. Similarly to the expression of proinflammatory cy- 
tokines, a significant increase in the IFN-B mRNA levels, and 
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FIG7 Characterization of rTGEV-mutG virus in cell culture. (A) Genetic structure of r(TGEV-mutG virus. The bar on top represents the TGEV genome, in 
which the different genes (ORF la, ORF 1b, S, 3a, 3b, E, M, N, and 7 genes), the leader sequence (L), the 3’ UTR, and the GAIT-like RNA motif are indicated. The 
structures depicted below the genome show the sequence and predicted secondary structure of the GAIT-like RNA motif from the WT (rTGEV-WT) and mutant 
(rTGEV-mutG) viruses, where the mutations introduced in the mutant are indicated in red. (B) In vitro growth kinetics of the WT and mutant viruses. ST cells 
were infected at an MOI of 0.05 (left) and 3 (right) with r—ITGEV-WT (WT) and rTGEV-mutG (mutG) viruses, and at different times postinfection the viral titer 
was determined by plaque titration on ST cells. Error bars represent the standard deviations from three independent experiments. 


those of its downstream expressed MDAS5 and RIG-I genes, was 
observed in cells infected with the rTGEV-mutG virus compared 
to those in cells infected with the WT virus (Fig. 9). These results 
indicated that the viral GAIT-like RNA motif interfered with the 
host innate immune response. 

Signaling pathway affected by the viral GAIT-like RNA mo- 
tif. In nonimmune cells, the production of cytokines and IFN is 
mediated mainly by RIG-I and MDAS5 cytoplasmic sensors. These 
proteins recognize dsRNAs, activating the transcription factors 
IRF3/7 and NF-«B, and the subsequent expression of proinflam- 
matory cytokines and IFN (16, 34). In the case of CoVs, it has been 
described that the viral RNA is recognized by MDAS (18), al- 


TGEV F2.2L 
kDa | WT mutG - 
170— Soe a EPRS 
2 — RRS 


FIG 8 Analysis of the interaction of EPRS and RRS with the mutated TGEV 
GAIT-like RNA motif. Proteins (600 wg) from cytoplasmic extracts of TGEV- 
infected Huh-7 cells were pulled down with either the native (F2.2L-WT) or 
the mutated (F2.2L-mutG) TGEV GAIT-like RNA motif, together without 
RNA (—) as the specificity control. The presence of EPRS and RRS in the initial 
extract (I) and in the pulled-down fractions was further analyzed by Western 
blotting. The molecular masses in kDa of EPRS and RRS are indicated at the 
left. 
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though a role for RIG-I was also proposed in several cell types (35). 
To analyze whether MDAS or RIG-I signaling was affected by the 
GAIT-like RNA motif, both sensors were independently silenced 
in ST cells by using specific small interfering RNAs (siRNAs). The 
silenced cells were then infected with r[GEV-WT or rTGEV- 
mutG viruses, and the effect on TNF, IFN-f, and viral genomic 
RNA production was analyzed by RT-qPCR in comparison with 
the levels in cells transfected with a negative-control siRNA. Two 
siRNA transfection steps were needed to achieve a sustained si- 
lencing at the mRNA and protein levels of about 43% and 35% for 
RIG-I and MDAS, respectively (Fig. 10A). Despite the modest 
reduction of MDAS expression, it was previously described that a 
similar reduction is enough to observe an effect on cytokine pro- 
duction (36). In infected MDA5-silenced cells but not in RIG-I- 
silenced cells, an increment in the viral genomic RNA levels of the 
rTGEV-WT and rTGEV-mutG viruses was observed compared 
with nonsilenced cells, being this increment higher in the case of 
the mutant virus (Fig. 10B, left). As previously observed (see 
above), in nonsilenced cells, infection by the rTGEV-mutG virus 
led to increased TNF and IFN-B mRNA levels compared with 
those caused by infection by the WT virus (Fig. 10B, right). Inter- 
estingly, in MDAS5-silenced cells, the rTGEV-mutG virus showed 
the same pattern of cytokine reduction as the WT virus, while this 
did not occur in RIG-I-silenced cells (Fig. 10B, right). These re- 
sults suggested that the GAIT-like RNA motif affects the MDA5 
signaling pathway. 
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FIG9 Accumulation kinetics of host genes involved in inflammation. ST cells were infected with either rTGEV-WT (blue) or rTGEV-mutG (red) viruses, and 
the expression of six genes involved in inflammation and innate immunity (TNF, CCL2, CXCL9, IFN-B, MDAS, and RIG-I) was analyzed at different hpi by 
RT-qPCR. In parallel, the expression of the TGF-B gene, whose expression is not affected by TGEV infection, and the levels of viral genomic RNA and sgmRNA 
7 were also evaluated as controls. The B-glucuronidase (GUSB) mRNA levels were used as an endogenous mRNA control in all cases. The mRNA levels were 
relative to those in noninfected cells. The fold change of viral RNA and cellular mRNA levels in cells infected with rTGEV-mutG compared with cells infected with 
rTGEV-WT is also indicated. Error bars indicate the standard deviations from three independent experiments. ***, P < 0.001; **, P< 0.01; *, P< 0.05. 


Inacomplementary approach, ST cells were treated with swine 
IFN-B for 16 h to induce the overexpression of RIG-I and MDA5 
(Fig. 11, left). Cells were then infected with the rTGEV-WT and 
rTGEV-mutG viruses, and the genomic RNA was quantified at 
8 hpi by RT-qPCR. A significant reduction in the genomic RNA 
levels of the rTGEV-WT virus was observed in cells treated with 
IFN-8 compared with those in untreated cells. Interestingly, this 
reduction was higher in IFN-f-treated cells infected with the 
rTGEV-mutG virus (Fig. 11, right) than in cells infected with the 
rTGEV-WT virus, confirming that the GAIT-like RNA motif 
modulates cytoplasmic sensor signaling. 


DISCUSSION 


The identification of host factors and the elucidation of their func- 
tions in viral replication are major areas in current virology re- 
search that will allow the development of novel antiviral strategies 
to control viral diseases. In this sense, we previously identified 
nine cellular proteins that interact with the 3’ end of the TGEV 
genome, from which a functional role in virus replication was 
described for hnRNP Q, PABP, and the aminoacyl-tRNA synthe- 
tases EPRS and RRS (13) (Marquez-Jurado et al., unpublished). In 
this study, we have identified a novel viral RNA motif, consisting 
of a 32-nt domain located 410 nt from the 3’ end of the TGEV 
genome, which directly interacted with the cellular proteins EPRS 
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and RRS. This RNA motif shared high homology in sequence and 
secondary structure with the cellular GAIT element that is in- 
volved in the modulation of the inflammatory response (29). Sim- 
ilar to the cellular GAIT element, we have shown that the viral 
RNA motif inhibited the in vitro translation of a chimeric mRNA 
containing this motif, indicating that in vitro it behaves as a GAIT- 
like RNA motif. Finally, a novel finding of this study was that the 
GAIT-like RNA motif counteracted the host innate immune re- 
sponse triggered by MDAS, at least in cell cultures. Further studies 
will be needed to confirm this effect in vivo and to analyze the role 
of the GAIT-like RNA motif in pathogenesis. 

In addition to the canonical function of aminoacyl-tRNA syn- 
thetases in translation, higher eukaryotic aminoacyl-tRNA syn- 
thetases have acquired a variety of noncanonical roles, including 
signaling, stress response, and transcriptional regulation (37, 38). 
Based in part on this versatility, these enzymes have been recruited 
by viruses to perform essential functions (39, 40). In this respect, 
EPRS and RRS could be involved in different essential processes of 
the TGEV life cycle, through their interaction with the 3’ end of 
the genome. 

A key observation of this study was the high homology in se- 
quences and secondary structures between the viral motif inter- 
acting with EPRS and RRS and the cellular GAIT element present 
at the 3’ UTR of several mRNAs encoding proinflammatory pro- 
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FIG 10 Effect of silencing the expression of RIG-I-like receptors on TNF and IFN-B production during rTGEV-WT and rTGEV-mutG infection. ST cells were 
transfected twice with a negative-control siRNA, an siRNA previously validated for porcine RIG-I, or an siRNA designed for porcine MDAS. After 24 h of the 
second transfection, cells were noninfected or infected with rTGEV-WT (WT) and rTGEV-mutG (mutG) viruses. Total RNA and proteins were extracted at 12 
and 18 hpi (36 and 42 h after the second siRNA transfection, respectively) to evaluate the silencing, virus production, and mRNA levels of TNF and IFN-B. (A) 
Silencing of RIG-I and MDAS. Silencing at the mRNA level (left) was evaluated at 36 (blue) and 42 (red) h after the second transfection of the siRNAs by RT-qPCR 
and expressed as the percentage of RIG-I and MDA5 mRNA remaining after the silencing of target genes in comparison to reference levels from cells transfected 
with the negative-control siRNA. RIG-I and MDAS silencing at the protein level (right) was evaluated by immunoblotting at 42 h after the second transfection 
of the RIG-I and MDAS siRNAs (siRIG-I and siMDAS, respectively) in comparison to cells transfected with the negative-control siRNA (siNEG). The B-actin was 
used as a loading control. (B) Effect of RIG-I and MDAS silencing on virus genomic RNA, IFN-B, and TNF production. Virus genomic RNA (gRNA) and IFN-B 
and TNF mRNA levels were quantified by RT-qPCR in cells transfected with the negative-control (siNEG; blue bars), RIG-I (siRIG-I; red bars), or MDA5 
(siMDAS5; green bars) siRNAs and infected with rTGEV-WT (WT) or rTGEV-mutG (mutG) viruses. Levels of virus gRNA are expressed as percentages in 
comparison with the gRNA levels in nonsilenced and infected cells (left). IFN-8 and TNF mRNA levels are relative to those in nonsilenced and noninfected cells 
(right). The figure shows the effects observed at 18 hpi (42 h after the second siRNA transfection), when they were more evident. In all cases, GUSB mRNA levels 
were used as the endogenous control. Errors bars indicate standard deviations from three independent experiments. ***, P < 0.001. 


teins (29). Monocytes and macrophages promote the synthesis of 
inflammatory cytokines and chemokines to counteract invading 
microorganisms. However, an unregulated synthesis of these in- 
flammatory products can have detrimental effects in the host. 
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FIG 11 IFN-B sensitivity of rTGEV-WT and rTGEV-mutG viruses. Conflu- 
ent ST cells were untreated (blue) or treated (red) with 4 ng/10° cells of porcine 
IFN-B (Kingfisher Biotech) for 16 h and then noninfected or infected with 
either rTGEV-WT (WT) or rTGEV-mutG (mutG) at an MOI of 3. At 8 hpi, 
total RNA was extracted, and the levels of viral genomic RNA (gRNA) were 
quantified by RT-qPCR and expressed relative to those in IFN-B-untreated 
and infected cells (right). As a control of the IFN- treatment, the MDA5 and 
RIG-I mRNA levels were quantified by RT-qPCR in noninfected cells, un- 
treated or treated with IFN-B (left). GUSB mRNA levels were used as an en- 
dogenous control in all cases. Error bars indicate the standard deviations from 
three independent experiments. *, P < 0.05. 
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Therefore, the host must have different mechanisms to regulate 
cytokine production and allow the resolution of inflammation. 
One of these mechanisms is provided by the GAIT system in 
monocytic cells (29). In response to IFN-y, the GAIT complex 
(constituted by EPRS, hnRNP Q, L13a, and GAPDH) is assembled 
and, by interaction with the GAIT element present in the 3’ UTR 
of several mRNAs encoding proinflammatory proteins, represses 
their translation. In this study, we have shown that the viral RNA 
motif is a GAIT-like RNA motif, based on its capacity to bind the 
GAIT complex and to inhibit the in vitro translation of a chimeric 
mRNA containing this motif. In this sense, EPRS and/or RRS 
could be components of a protein complex that regulate viral 
translation by its interaction with the viral GAIT-like RNA motif. 
Nevertheless, the characterization of an rTGEV-mutG virus, har- 
boring mutations that disrupt the GAIT-like RNA motif, indi- 
cated that the GAIT-like RNA motif was not critical in CoV trans- 
lation and replication in cell cultures, suggesting that during 
infection the viral RNA motif could have a function different to 
that described for the cellular GAIT element. 

After confirming the dispensability of the GAIT-like motif for 
virus replication in cell cultures, its potential role modulating the 
host inflammatory response and more specifically the innate im- 
mune response was further analyzed in cell cultures. The innate 
immune response is the first line of antiviral defense that triggers 
the synthesis of proinflammatory cytokines and type I IFN to pre- 
vent viral infections (15). In this study, an exacerbated expression 
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TABLE 1 Oligonucleotides used for PCR amplifications 


Oligonucleotide? 


Oligonucleotide sequence (5'—>3")?’ 


Xmal-sp6-5’-VS 
5'TGEV-504-RS 
Xmal-T7-VS 
TGEV-3'-RS 
TGEV-3'-pA-RS 


CGCGCCCGGGATTTAGGTGACACTATA(G)ACTTTTAAAGTAAAGTGAGTGTAGC 
CACCAATGACGTAGTGATCCTTACG 
CGCGCCCGGGTAATACGACTCACTATA(G)GGATGCTGTATTTATTACAGTTTTAATC 
GTATCACTATCAAAAGGAAAATTTTCAAACAATTTTAC 
TITTTTTTTTITTTTTTTTTTTTTTGTATCACTATCAAAAGG 


F1-RS 
F2-VS 
F2-RS 
F3-VS 
F3-RS 
F4-VS 
F4-RS 
F5-VS 
F5-RS 
F6-VS 


4 VS, virus sense; RS, reverse sense. 


CTATATAAGATATTAAAGTCATTGACAGTTTTAAGATTG 
CGCGCCCGGGTAATACGACTCACTATA(G)GGACTATTACTTAATCACTCTTTCAAT 
GCAAAATCCTAGTAAGACTAGAAAGATTACTCGAAGC 
CGCGCCCGGGTAATACGACTCACTATA(G)GGATTTTGCTGCTACAGATTGTTAG 
GTGCTTACCATTCTGTACAAGAGTAGACAGCG 
CGCGCCCGGGTAATACGACTCACTATA(G)GGCTGTCTACTCTTGTACAGAATGGTA 
CTTTAAACTTCTCTAAATTACTAAATCTAGCATTGC 
CGCGCCCGGGTAATACGACTCACTATA(G)GGTTAGTCACATTAATGTAAGGCAAC 
CTAATATGCAATAGGGTTGCTTGTACCTCCTATTAC 
CGCGCCCGGGTAATACGACTCACTATA(G) GGAGAAGTTTAAAGATCCGCTACGA 


’ Transcription promoters are in bold. Transcription initiation sites are in parentheses. The template sequence corresponding to poly(A) is underlined 


of both proinflammatory cytokines and IFN-f, together with its 
downstream ISGs (MDA5 and RIG-I), has been observed during 
infection with the rTGEV-mutG virus, harboring mutations in the 
GAIT-like RNA motif, compared with that in the rTGEV-WT 
virus infection. The level of the effect observed on the expression 
of proinflammatory cytokines and IFN-6 was similar or even 
higher than that reported for CoV and other RNA viruses (14). In 
addition, we have identified that the MDA5 pathway was the in- 
nate immune pathway significantly affected by the GAIT-like 
RNA motif. In contrast to other RNA viruses, in which structured 
RNA domains present in their genome are recognized by the cell 
sensor triggering the innate immune response (20-23, 41), the 
data presented in this study suggested that the TGEV GAIT-like 
RNA motif interferes with the innate immune response, at least in 
cell cultures. Regarding the mechanism of action of the GAIT-like 
RNA motif, it could be postulated that this motif may act directly, 
preventing the activation of MDAS. Alternatively, the interaction 
with EPRS and RRS could mask the GAIT-like RNA motif to make 
it inaccessible to cellular sensors, such as MDAS, thus avoiding the 
triggering of the innate immune response. In this sense, mutations 
in the viral GAIT-like RNA motif could facilitate its sensing by 
MDAS, eliciting a faster and stronger antiviral response, such as 
the one observed during the infection with the rTGEV-mutG vi- 
rus. Similar mechanisms have been described in negative- 
stranded RNA viruses, such as influenza and Ebola viruses, in 
which a viral protein instead of a cellular protein binds a dsRNA 
signature in their genomes, hiding them from recognition by the 
cellular sensors (14). To our knowledge, this is the first report ofa 
cellular protein proposed to mask a viral PAMP RNA domain to 
subvert the host antiviral response. 

The identification of the TGEV GAIT-like motif and the clar- 
ification of its role in the modulation of the antiviral response may 
be of relevance for other coronaviruses, such as the highly patho- 
genic SARS-CoV and MERS-CoV. In fact, several GAIT-like RNA 
motif candidates have been found within the 3’ UTR of SARS- 
CoV and MERS-CoV by bioinformatics analysis (unpublished re- 
sults), although complementary approaches will be needed to 
clarify whether these putative GAIT-like RNA motifs could act as 
such. 
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MATERIALS AND METHODS 


Cells and viruses. Swine testis (ST) cells (ATCC, CRL-1746) (42) were 
grown in Dulbecco’s modified Eagle medium (DMEM) supplemented 
with 10% fetal bovine serum (FBS). Baby hamster kidney cells, stably 
transformed with the porcine aminopeptidase N gene (BHK-pAPN; 
kindly provided by H. Laude, Unité de Virologie Immunologie Molécu- 
laires, France) to make them susceptible to TGEV, were grown in DMEM 
supplemented with 5% FBS and Geneticin (1.5 mg/ml) as the selection 
agent. Human liver-derived Huh-7 cells (kindly provided by R. Barten- 
schlager, University of Heidelberg, Germany) were grown in DMEM sup- 
plemented with 10% heat-inactivated FBS. Human monocytic U937 cells 
(kindly provided by M. Fresno, Centro de Biologia Molecular Severo 
Ochoa, CBMSO, CSIC-UAM, Spain) were grown in RPMI 1640 medium 
supplemented with 10% FBS. 

The TGEV PUR46-MAD strain (43) was used to infect ST cells, and 
the TGEV PUR46-C11 strain (44) was used to infect Huh-7 cells. In all 
cases, virus titration was performed on ST cell monolayers as previously 
described (45). 

DNA constructs and RNA transcripts used for RNA affinity chroma- 
tography. The DNA templates for in vitro transcription of the precise 
TGEV 5’ end, TGEV 3’ end without or with a poly(A) tail, and fragments 
Fl to F6 were generated by PCR by using specific oligonucleotides (Ta- 
ble 1) and the pML33L-637G plasmid (46) as the template. To amplify the 
first 504 nt of the TGEV 5’ end, PCR was performed with the oligonucle- 
otides XmalI-SP6-5’-VS, which includes the SP6 promoter, and 5’'TGEV- 
504-RS. To amplify the last 493 nt of the genome without or with a 25-nt 
poly(A) tail, the PCRs were performed with the oligonucleotides Xmal- 
T7-VS (including the T7 promoter) and either TGEV-3’-RS to generate 
the 3’ end without poly(A) or TGEV-3’-pA-RS to amplify the 3’ end with 
the poly(A) tail. Finally, to amplify the 3’-end fragments F1 to F6, the 
PCRs were performed with the oligonucleotides XmaI-T7-VS and F1-RS 
for fragment F1, F2-VS and F2-RS for fragment F2, F3-VS and F3-RS for 
fragment F3, F4-VS and F4-RS for fragment F4, F5-VS and F5-RS for 
fragment F5, and F6-VS and TGEV-3’-RS for fragment F6. In all cases, the 
virus sense primer (VS) included the T7 promoter to allow the in vitro 
transcription. All PCRs were performed with the platinum Pfx DNA poly- 
merase (Invitrogen) by following the manufacturer’s recommended con- 
ditions. The resulting amplicons were purified with the QIAquick gel 
extraction kit (Qiagen) and used as the templates for in vitro transcription. 

In vitro transcription reactions were performed from 600 ng of each 
PCR amplicon using the SP6/T7 transcription kit (Roche). To label the 
transcript with biotin, biotin-14-CTP (Invitrogen) was added to the tran- 
scription reaction at a final concentration of 0.16 mM in a 1:6.25 ratio to 
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TABLE 2 Biotin-labeled RNAs used in affinity chromatography experiments 


RNA RNA bait sequence (5’—3’) 

F2.1 AGACUAUUACUUAAUCACUCUUUCAUACUAAAAGUGGUGCUUCGAGUAAU. 

F2.2 AUCUUAAAACUGUCAAUGACUUUAAUAUCUUAUAUAGGAGUUUAGCAGAAACCAGAU 
F2.3 CUUUCUAGUCUUACUAGGAUUUUGCU 

F2.2U GACUUUAAUAUCUUAUAUAGGAGUU 

F2.2D AAUCUUAAAACUGUCAAUGUUAGCAGAAACCAGAUU 

F2.2L AAUCUUAAAACUGUCAAUGACUUUAAUAUCUU 

F2.2L-mutG AAUCUUAAGACCGUCAAUGAUUUUAAUAUCUU 

F2.2R AUAUAGGAGUUUAGCAGAAACCAGAUU. 


unlabeled CTP. The transcription reaction mixtures were incubated for 
2 h at 37°C and treated with 10 U of DNase I (Roche) for 30 min at 37°C. 
The resulting transcripts were purified with the RNeasy minikit (Qiagen) 
by following the RNA cleanup protocol, analyzed by denaturing electro- 
phoresis in 2.5% (wt/vol) agarose gels, and spectrophotometrically quan- 
tified. 

The RNA fragments F2.1, F2.2, F2.3, F2.2U, F2.2D, F2.2L, F2.2L- 
mutG, and F2.2R (Table 2), labeled at the 3’ end with biotin, were chem- 
ically synthesized (Sigma Aldrich). 

RNA affinity chromatography. Huh-7 cells were grown to confluence 
in 15-cm-diameter dishes and infected at an MOI of 5 with TGEV PUR46- 
C11. At 16 hpi, cytoplasmic cell extracts were prepared as described pre- 
viously (13), and total protein concentration was determined with the 
BCA protein assay kit (Thermo Scientific). To perform the RNA affinity 
chromatography, 600 yg of cytoplasmic cell extract was diluted 1:1 in 
binding buffer (50 mM HEPES [pH 7.9], 150 mM KCl, 5% glycerol, 
0.01% NP-40) and precleared three times by incubation with 20 pl of 
streptavidin Sepharose beads (streptavidin Sepharose high performance; 
Amersham) for 4 h at 4°C. Biotinylated RNA (between 5 and 10 wg) was 
diluted in 200 yl of RNA-binding buffer (5 mM HEPES [pH 7.9], 0.5 mM 
EDTA, 1 M KCl) and incubated with 20 pl of fresh streptavidin Sepharose 
beads at 4°C for 1 h. The immobilized RNA was washed three times for 
15 min with 200 pl of binding buffer and then incubated with the pre- 
cleared protein extract overnight at 4°C. The RNA-protein complexes 
were washed three times for 15 min with 200 pl of binding buffer, pulled 
down, and resuspended in 30 jl of NuPAGE sample buffer (Invitrogen). 
After incubation for 10 min at 65°C, the beads were centrifuged and the 
supernatant containing RNA-interacting proteins was recovered. All in- 
cubations were performed at 4°C in an orbital shaker, and the beads were 
collected by centrifugation at 500 X g. 

Western blot analysis. Cell lysates or pulled down RNA-binding pro- 
teins were resolved by denaturing electrophoresis in NuPAGE 4 to 12% 
bis-Tris gels with morpholinepropanesulfonic acid (MOPS) SDS running 
buffer (Invitrogen) and transferred electrophoretically to nitrocellulose 
membranes (Hybond-C extra nitrocellulose; Amersham Biosciences) by 
following the manufacturer’s recommended conditions. Membranes 
were blocked for 1 h at room temperature with 5% dried skimmed milk in 
TBS (20 mM Tris-HCl [pH 7.5], 150 mM NaCl) and then probed with 
specific antibodies diluted in TBS and 0.1% Tween 20 (TTBS) containing 
3% dried skimmed milk. Primary antibodies for hnRNP A1 (rabbit poly- 
clonal E17), hnRNP Q (mouse monoclonal antibody [MAb] I8E4:sc- 
56703), and hnRNP AO (goat polyclonal G-17:sc-16509) were purchased 
from Santa Cruz Biotechnology. Antibodies for hnRNP U (rabbit poly- 
clonal ab20666), EPRS (rabbit polyclonal ab31531), RRS (rabbit poly- 
clonal ab31537), and B-actin (mouse MAb ab8226) were purchased from 
Abcam. The p100 transcriptional coactivator guinea pig polyclonal anti- 
serum GP25 was purchased from Progen Biotechnik. The MDAS5 rabbit 
monoclonal antibody (D74E4) and the RIG-I rabbit polyclonal antiserum 
were purchased from Cell Signaling and AnaSpec, respectively. The PABP 
rabbit polyclonal antiserum was kindly provided by Amelia Nieto (Centro 
Nacional de Biotecnologia, CNB-CSIC, Spain). The blots were then incu- 
bated with horseradish peroxidase-conjugated secondary antibodies di- 
luted in TTBS-3% dried skimmed milk for 1 h at room temperature, and 
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the immune complexes were detected using the Immobilon Western 
chemiluminescent substrate (Millipore) according to the manufacturer’s 
instructions. 

RNA EMSA. RNA-protein binding reactions were performed by incu- 
bating 100 pmol of the 3’ biotinylated TGEV GAIT-like RNA motif (Ta- 
ble 2, F2.2L) with 1 pmol of recombinant purified EPRS (OriGene), 
GAPDH (Abcam), and DDK/Myc tag (NH2-EQKLISEEDLAANDILDY 
KDDDDKV-OH) in binding buffer (12% glycerol, 20 mM Tris-HCl 
[pH 7.4], 50 mM KCI, 1 mM EDTA, 1 mM MgCl, 1 mM dithiothreitol 
[DTT]) containing 0.2 U/l SUPERase inhibitor (Ambion) for 30 min at 
25°C. The RNA-protein complexes were resolved in a 4% nondenaturing 
PAGE gel using 0.5 TBE as the running buffer. After electrophoresis, the 
gel was blotted onto positively charged nylon membranes (BrightStar- 
Plus; Ambion), and the biotinylated RNA was detected using the Bright- 
Star BioDetect kit (Ambion), following the manufacturer’s instructions. 

RNA immunoprecipitation assay. ST cells were grown in 15-cm- 
diameter dishes to confluence and infected at an MOI of 5 with TGEV 
PUR46-MAD. At 8 hpi, the cells were harvested and centrifuged at 
2,000 X g for 5 min at 4°C, and the cell pellets were resuspended in RIP 
lysis buffer (MagnaRIP; Millipore) containing protease inhibitor cocktail 
(Roche) and 1 U/pl RNasin RNase inhibitor (Promega). Then, the cell 
suspension was incubated at 4°C for 5 min and centrifuged at 3,000 X g for 
2 min at 4°C, and the supernatant was collected to perform the RNA 
immunoprecipitation assay using the EZ-Magna RIP kit (Millipore) by 
following the manufacturer’s instructions. Briefly, cytoplasmic superna- 
tant, corresponding to 1 X 107 cells, was immunoprecipitated with anti- 
EPRS, anti-RRS, or anti-GFP (Roche) antibodies (attached to magnetic 
beads). Immunoprecipitated RNA-protein complexes were incubated 
with proteinase K for 30 min at 55°C, and the RNA was isolated using the 
RNeasy minikit (Qiagen), according to the manufacturer’s recommenda- 
tions. Finally, the purified RNA was subjected to RT-qPCR for the detec- 
tion of EPRS-, RRS-, or GFP-associated viral genomic RNA and sgmRNAs 
by using specific oligonucleotides (Table 3) and SYBR green-based detec- 
tion. 

Generation of plasmids and RNA constructs for in vitro translation. 
To analyze the effect on translation of the TGEV GAIT-like RNA motif, 
the plasmids pINT-Luc-Cp, pTNT-Luc-TGEV, and pTNT-Luc were 
generated. These plasmids expressed (under the control of the T7 pro- 
moter) three chimeric mRNAs carrying the luciferase gene, which con- 
tained at the 3’ UTR the Cp GAIT element, the TGEV GAIT-like RNA 
motif, or no GAIT element. The pTNT-Luc-Cp plasmid was generated by 
cloning into the pTNT vector (Promega) digested with Mlul and NotI a 
chemically synthesized (GeneArt; Life Technologies) DNA fragment 
(1,761 bp), which contained the restriction site Mlul, the luciferase gene, 
the first 48 nt of the luciferase mRNA 3’ UTR, a 15-nt polylinker, and the 
29-nt Cp GAIT element (28), flanked by the restriction sites Sall and Notl. 
To generate the pT[NT-Luc-TGEV plasmid, the SalI-NotI DNA fragment 
from pINT-Luc-Cp was exchanged with a chemically synthesized (Ge- 
neArt; Life Technologies) SalI-Not] DNA fragment containing the 32-nt 
TGEV GAIT-like RNA motif. Finally, the plasmid pTNT-Luc, which did 
not contain any GAIT element, was generated by religation of pTNT- 
Luc-Cp digested with SalI and NotI and treated with Klenow polymerase. 

To generate the capped mRNA construct for in vitro translation, the 
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TABLE 3 Oligonucleotides used for RT-qPCR of viral gRNA and sgmRNA 


Amplicon Oligonucleotide* Oligonucleotide sequence (5'>3")’ Method 
gRNA Rep-VS TTCTTTTGACAAAACATACGGTGAA SYBR 
Rep-RS CTAGGCAACTGGTTTGTAACATCTTT SYBR 
sgmRNA S LCS1-VS CCAACTCGAACTAAACTTTGGTAACC SYBR 
LCS1-RS TCAATGGCATTACGACCAAAAC SYBR 
sgmRNA 3a rt3a-VS2 CGGACACCAACTCGAACTAAACTTAC SYBR 
rt3a-RS ATCAAGTTCGTCAAGTACAGCATCTAC SYBR 
sgmRNA E Ldrt-VS CGTGGCTATATCTCTTCTTTTACTTTAACTAG SYBR 
mRNAE-RS CAATGCCCTAGGAAACGTCATAG SYBR 
sgmRNA M Ldrt-VS CGTGGCTATATCTCTTCTTTTACTTTAACTAG SYBR 
mRNAM-RS GCATGCAATCACACACGCTAA SYBR 
sgmRNA N Ldrt-VS CGTGGCTATATCTCTTCTTTTACTTTAACTAG SYBR 
N82—RS TCTTCCGACCACGGGAATT SYBR 
sgmRNA 7 Ldrt-VS CGTGGCTATATCTCTTCTTTTACTTTAACTAG SYBR 
7 (38)-RS AAAACTGTAATAAATACAGCATGGAGGAA SYBR 
gRNA Rep-VS TTCTTTTGACAAAACATACGGTGAA TaqMan 
Rep-RS CTAGGCAACTGGTTTGTAACATCTTT TaqMan 
Rep-MGB FAM-AGGGCACCGTTGTCA-MGB TaqMan 
sgmRNA 7 Ldrt-VS CGTGGCTATATCTCTTCTTTTACTTTAACTAG TaqMan 
7 (38)-RS AAAACTGTAATAAATACAGCATGGAGGAA TaqMan 
mRNA7-MGB FAM-CGAACTAAACGAGATGCT-MGB TaqMan 


“VS, virus sense; RS, reverse sense; MGB, minor groove binder group. 
» FAM, 6-carboxyfluorescein. 


© SYBR, SYBR green-based detection chemistry for qPCR; TaqMan, TaqMan-based detection for qPCR. 


pINT-Luc-Cp, pINT-Luc-TGEV, and pTNT-Luc plasmids were linear- 
ized with BamHI and in vitro transcribed with T7 polymerase by using the 
mMESSAGE mMACHINE 17 kit (Ambion), by following the manufac- 
turer’s recommendations. The resulting transcripts were purified with the 
RNeasy minikit (Qiagen), analyzed by denaturing electrophoresis in 2.5% 
(wt/vol) agarose gels, and spectrophotometrically quantified. 

Preparation of cytosolic extracts from U937 cells. To prepare cyto- 
solic extracts of U937 cells, 1 X 107 cells grown in 10-cm-diameter dishes 
were incubated in 10 ml of RPMI 1640 containing 0.5% FBS for 1 h at 
37°C and then treated or not treated with 500 U/ml of recombinant hu- 
man IFN-y (R&D Systems) for 24 h at 37°C. After IFN treatment, cyto- 
solic extracts were prepared using the PhosphoSafe extraction reagent 
(Novagen) containing protease inhibitor cocktail (Roche), by following 
the manufacturer’s recommendations, and the protein concentration was 
determined using the bicinchoninic acid (BCA) protein assay kit (Thermo 
Scientific). 

In vitro translation of chimeric luciferase mRNAs. In vitro tran- 
scribed chimeric luciferase mRNAs, containing the Cp GAIT element, the 
TGEV GAIT-like RNA motif, or no GAIT element at the 3’ UTR, were in 
vitro translated in the absence or in the presence of the GAIT complex. To 
this end, 200 ng of capped chimeric transcripts were in vitro translated 
with 35 yl of rabbit reticulocyte lysate (Promega) in the presence of 4 wg 
of cytosolic extracts from U937 monocytic cells untreated or treated with 
IFN-y, 1 ~M amino acid mixture, and 40 U of RNasin RNase inhibitor 
(Promega). The translation reactions were performed at 30°C for 90 min 
in a final volume of 50 yl. Finally, the translation efficiency was analyzed 
by luciferase activity quantification using the luciferase assay system 
(Promega) and an Orion microplate luminometer (Berthold Detec- 
tion System). 

Construction of a TGEV cDNA clone harboring mutations in the 
GAIT-like RNA motif. The TGEV infectious clone harboring mutations 
in the GAIT-like RNA motif (pBAC-TGEV-mutG) was generated from a 
TGEV cDNA clone (pBAC-TGEV?!-S,,,,) containing the S gene of the 
TGEV Purdue-type virus isolate (32, 44, 47). To this end, a DNA fragment 
(mutG) expanding the last 727 nt of the viral genome was generated by 
chemical synthesis (GeneArt; Life Technologies). This fragment was 
flanked by the restriction sites KpnI and BamHI and contained three point 
mutations, A to G, U to C, and C to U, at genomic positions 28178, 28181, 
and 28190, respectively. The DNA fragment was digested with KpnI and 
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BamHI and cloned into the intermediate plasmid pSL-3EMN7 (contain- 
ing the TGEV genome nt 20372 to 28087) digested with the same restric- 
tion sites, leading to plasmid pSL-3EMN7-mutG. To generate pBAC- 
TGEV-mutG, the Mlul-BamHI DNA fragment from the infectious clone 
pBAC-TGEVF"-S,y was exchanged with the Mlul-BamHI DNA frag- 
ment from plasmid pSL-3EMN7-mutG. All cloning steps were checked by 
restriction analysis and sequencing. 

Recovery of the recombinant virus from the cDNA clone. To recover 
the rTGEV-mutG recombinant virus, BHK-pAPN cells were grown to 
95% confluence on 35-mm_-diameter plates and transfected with 4 wg of 
infectious cDNA clone pBAC-TGEV-mutG using 12 pl of Lipofectamine 
2000 (Invitrogen), according to the manufacturer’s specifications. After 
6h of incubation at 37°C, cells were trypsinized, plated over confluent ST 
monolayers grown on a 35-mm-diameter plate, and incubated at 37°C for 
72 h. The cell supernatants were harvested, and the recovered virus was 
cloned by three rounds of plaque purification. Recombinant virus was 
grown and titrated as previously described (45). 

Analysis of cellular gene expression and viral RNA synthesis. Cellu- 
lar gene expression and viral RNA synthesis were quantified by RT-qPCR. 
Total RNA from uninfected or infected ST cells was prepared using the 
RNeasy minikit (Qiagen), and total cDNA was synthesized with random 
hexamers from 100 ng of total RNA, using the high-capacity cDNA 
transcription kit (Life Technologies), by following the manufacturers’ 
instructions. Cellular gene expression was analyzed using TaqMan 
gene expression assays (Life Technologies) specific for porcine TNF 
(Ss03391318_g1), CCL2 (Ss03394377_m1), CCL4 (Ss03381395_ul), 
CXCL9 (Ss03390033_m1), CXCL11 (Ss03648935_g1), IFN-B 
(Ss03378485_ul), MDA5 (Ss03386373_ul1), RIG-I (Ss03381552_ul), 
or TFG-B (Ss03382325_u1) genes. To analyze viral RNA levels, custom 
TaqMan assays (Life Technologies) specific for the TGEV genomic 
RNA and sgmRNA 7 (Table 3) were used. The B-glucuronidase 
(GUSB) gene (TaqMan code $s03387751_u1) was used as a reference 
housekeeping gene, since its expression remains constant in infected 
cells compared to that in uninfected cells. Data were acquired with a 
7500 real-time PCR system (Applied Biosystems) and analyzed with 
7500 software v2.0.6. All experiments and data analysis were MIQE 
compliant (48). 

siRNA transfection and silencing analysis. ST cells were transfected 
twice with specific Silencer Select siRNAs (Ambion-Life Technologies) by 
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following a reverse transfection protocol. Briefly, for each well of a 6-well 
plate, 8.4 X 10° cells were transfected in suspension with 25 nM of each 
siRNA using 8 pl of siPORT amine (Ambion-Life Technologies) diluted 
in 100 wl of Opti-MEM medium (Gibco BRL-Life Technologies), by fol- 
lowing the manufacturer’s instructions. For the porcine MDAS5 gene, the 
siRNA used in the assay (sense sequence of 5’ CCC AGA AUA UGA AAA 
UGA Att 3’ and antisense sequence of 5’ UUC AUU UUC AUA UUC 
UGG Gtt 3’ [lowercase indicates nucleotides protruding at the 3’ ends]) 
was selected from three designed siRNAs, as it provided the highest silenc- 
ing efficiency (around 40%) in previous MDA5-silencing experiments 
(data not shown). To silence the RIG-I gene, an siRNA (sense sequence of 
5' CAG UCA GAG ACA ACU UGG Att 3’ and antisense sequence of 
5’ UCC AAG UUG UCU CUG ACU Gtt 3’ [lowercase indicates nucleo- 
tides protruding at the 3’ ends]) previously validated in porcine cells was 
used (36). As a negative control, an irrelevant siRNA (sense sequence of 
5’ CUGCCCCAGCGAUAUCCAGtt 3’ and antisense sequence of 5’ CU 
GGAUAUCGCUGGGGCAGtt 3’ [lowercase indicates nucleotides pro- 
truding at the 3’ ends]) previously described for porcine cells (36) was 
used. After transfection, cells were plated onto each well using DMEM 
with 10% FBS without antibiotic and incubated at 37°C for 48 h. Then, 
cells were trypsinized, retransfected in suspension with 25 nM of each 
siRNA, seeded in 24-well plates at a confluence of 2 X 10° cells per well, 
and incubated at 37°C for 24 h. Transfection medium was then discarded, 
and cells were infected with rTGEV-WT and rf[GEV-mutG viruses at an 
MOI of 3. Total RNA and protein extract were collected at 12 and 18 hpi 
(36 and 42 h after the second siRNA transfection, respectively) for further 
analysis. 

Statistical analysis. A two-tailed, unpaired Student f test was used to 
analyze differences in mean values between groups. All results were ex- 
pressed as means + standard deviations of the means. P values of <0.05 
were considered significant. 
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